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Using a field-theoretical representation of the Tanaka-Edwards integral [J we develop a method to 
systematically compute the number Na of 1-spin-stable states (local energy minima) of a glassy Ising 
system with nearest-neighbor interactions and random Gaussian couplings on an arbitrary graph. 
In particular, we use this method to determine Ng for K-vegulax random graphs and d-dimensional 
regular lattices for d — 2, 3. The method works also for other graphs. Excellent accuracy of 
the results allows us to observe that the number of local energy minima depends mainly on local 
properties of the graph on which the spin glass is defined. 



I. INTRODUCTION 

Glassy systems have non-trivial energy landscapes just 
as many complex systems observed in nature. The main 
characteristics of such landscapes is the number Ns of 
local minima, called also metastable states. Typically, 
this number grows exponentially with the system size N: 
Ns ~ e^f'. The rate of the exponential growth /* is 
a fundamental quantity characterizing the complexity of 
the system. It is however very difficult to calculate, and 
it has been analytically found in only a few cases: for the 
Ising model with random Gaussian interactions on a com- 
plete graph (SK-model) 0, on the one-dimensional 
closed chain [sl, f3j and for the Ising model with random 
binary interactions J = ±1 on if- regular random graphs 

i. 

In this paper we describe a method to determine /* 
for the Ising model with random Gaussian interactions 
on an arbitrary graph. In particular we use this method 
to compute /* for some JsT- regular graphs. The idea is 
to express the number of metastable states Ng in terms 
of the Tanaka-Edwards integral [ij and then to treat this 
expression as the partition function of a certain statisti- 
cal field theory. The logarithm of the partition function 
can be represented as a sum of connected Feynman dia- 
grams which in turn can be generated and summed on a 
computer, up to a certain order of the perturbative series. 
Using some general properties of this series we are able 
to estimate the value /, already from the first few orders 
with a very good accuracy. This allows us to observe 
that the values of /* are very similar for regular graphs 
with different topologies, indicating that the number of 
metastable states depends mainly on local properties of 
the graph. 



II. DERIVATION OF THE STATISTICAL FIELD 
THEORY 

We consider a system of Ising spins, cr, = ±1, i = 
1, . . . , N, residing on nodes of a simple graph described 
by an adjacency matrix A. The graph does not need to 
be connected. The matrix A is an x symmetric 



matrix with Aij — 1 if i and j are connected by an edge 
or Aij — otherwise. The energy of the system is given 
by: 
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where the coupling constants Jij's are random numbers 
taken from the standardized Gaussian distribution with 
zero mean and unit variance. 

We are interested in counting the number of local min- 
ima Ns of the energy ^ . A local (metastable) minimum 
is defined as a configuration of spins {ai} such that a fiip 
of any single spin increases energy. Such a configuration 
is also called 1-spin-stable. The number of 1-spin-stable 
states is given by [J, O] 

Ns= J2 ••• E I[(^Uj2^^^J^i''X^ (2) 
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where 6{x) is a step (Heaviside) function. In Ref. [l| it 
was shown that the averaging over Gaussian couplings 
Jij leads to the following concise formula: 

{Ns) = / Hl^qke-^^^-^^'^^'^^''^, (3) 
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Here ki = Ay is called degree of node i. The inte- 
gral ([3]) was calculated pj in the limit of oo for a 
complete graph: Aij = 1 — 5ij , yielding the result 



In {Ns)^Nf,, 



(6) 



with /* « 0.199228, known also from earlier considera- 
tions of the SK model 
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In this paper we shall propose a systematic method 
to evaluate this integral also for other graphs. Let us 
introduce an auxihary constant g to Eq. ([3]): 
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The idea is now to find a systematic way of expanding 
([7]) in powers of g and then to use this series expansion 
to estimate its value for 5 = 1: 



(7V.) = (iV,) (5)1,^1. 



(8) 



Borrowing some techniques from field theory let us define 
the following generating function: 



Z[J]= / \l^q^e'T..J.<i.^ 

k=i 



which allows for rewriting Eq. ([7| as: 



(iV,) (g) = exp 
The function Z[J] has a closed form: 
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where erf(a;) is the error function and c„ are cumu- 
lants of 1 + erf (jfc/\/2). The coefficients c„ can be 
easily calculated up to an arbitrary or der using a pro- 
gram for symbolic calculations: ci — ^/S/tt, C2 = — S/tt, 
C3 = 2(4 - 7r)/(%/27r3/2), ... . Equation ^ can be 
graphically represented as a sum of vacuum Feynman 
diagrams of a field theory with the propagator gMij and 
<i>"-vertices with coupling constants c„. The Feynman 
rules to calculate the contribution of a diagram are as 
follows. To each vertex at which n lines meet, we 
ascribe a factor c„. The subscript j means that the ver- 
tex is decorated by an index j = 1, . . . ,N which can be 
thought of as a color taken from a palette of N possible 
colors. A line joining two vertices decorated with colors 
i and j contributes a factor gMij . Additionally each di- 
agram has a certain symmetry factor which depends on 
the shape of the diagram. Finally, one needs to perform 
the summation over colors. 

As usual, the logarithm of Eq. (fTO|) . 



F{g) = \n{N,) (g), 



(12) 



contains only the contribution from connected diagrams. 
We shall see below that F{g) is an extensive quantity in 
N. Thus it is convenient to introduce a density Foo{g) 
per spin which for large N becomes a function of g only 
and can be represented as a power series: 
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Our goal is to determine /* = i^oo(l), which gives the 
rate of exponential growth of the number of 1-spin-stable 
configurations for large N. 

The coefficients at g^ in the series expansion (fT3|) come 
from connected Feynman diagrams with I links. In the 
general case, they must be summed on a computer, be- 
cause their number grows very fast with I. We need a 
systematic procedure which allows us to sum diagrams 
in order to calculate the coefficients /;. In our approach, 
such a procedure looks as follows: (A) "draw" all possible 
connected Feynman diagrams with I links; (B) calculate 
their symmetry factors s; (C) decorate each vertex of the 
diagram with an index i — 1, . . . , N; (D) for each deco- 
ration calculate the contribution of the diagram to F{g) 
as: 
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where v is the number of vertices of the diagram, the 
first product goes over all vertices and the second one 
over all links of the diagram; (E) sum contributions of all 
diagrams and decorations. 

Clearly, the procedure described above is not efficient 
if the original graph is sparse since then many propa- 
gators gMij are zero. Therefore, many of possible TV" 
decorations give a zero contribution and one wastes time 
summing many zeros. In particular, all decorations of a 
diagram having a self-connecting link give no contribu- 
tion since Ma — 0. Thus one can omit such diagrams in 
the sum. One can improve the step (C) of the procedure 
by concentrating only on decorations which potentially 
have a chance to contribute. In other words, one should 
look only for decorations for which propagators Mij do 
not vanish. This means that i and j must be neighbors on 
the graph on which the spins reside. Therefore, instead of 
the step (C) one should take the step (C) in which one 
only checks such decorations which are consistent with 
the graph structure. This can be done iteratively. First, 
one assigns a label j = 1, . . . , iV to one vertex of the di- 
agram. Then one assigns to its neighbors only values i 
such that Mij 7^ or equivalently that j and i are neigh- 
bors on the original graph. One repeats this process for 
neighbors of z's and so on, and selects only those labels 
for which the propagator does not vanish. This speeds 
up the step (E) of the procedure since now the number 
of decorations is of order Nk^~^ , where k is the average 
node degree of the graph. 

The procedure (A-C'-E) works for any graph. For a 
if-regular graph one can essentially simplify calculations 
since in this case Mij = Aij/K, as it stems from Eq. |[5|), 
and thus the elementary contribution (fT4l) to F{g) is 



(15) 



The last product of A^-'s is either zero or one. Summing 
over all decorations of a given diagram we get some num- 
ber P of decorations consistent with the graph structure. 
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FIG. 1: Feynman diagrams for I < 4, without multiple links, 
generated by Eqs. ifTOll and ifTTjl. together with their symme- 
try factors. 



Each Feynman diagram has now the following contribu- 
tion: 



P 



(16) 



Actually P is the only part of the expression which de- 
pends on the graph structure on which spins reside, other 
quantities can be calculated beforehand. Therefore, we 
used a C++ program to generate all simple diagrams 
(without multiple- and self-connections) up to ^ = 13 
and appHed the routine nauty IQ] for isomorphism test- 
ing to distinguish different diagrams and to determine 
their symmetry factors s. In figure [T] we show several ex- 
amples of diagrams for I < 4. Their number grows very 
fast with I: for I = 11 there are 11461 such diagrams, 
for I ^ 12 - 40964 and for ^ = 13 - 153786. The num- 
ber of all multi-diagrams, that is diagrams with multiple 
connections, is much larger but luckily there is no need 
to generate them. Any multi-diagram can be obtained 
from a simple diagram 7 by replacing its links by multi- 
ple links with a certain multiplicity m. The contribution 
of all multi-diagrams associated with a simple diagram 7 
is: 



00 

•E 

mi = l 



(17) 



where s and P are calculated for 7. The terms in the 
sum contribute to the order /' = mi + . . . + m; of the 
expansion of F{g), and V is the total number of links 
of the multi-diagram. Factors l/nial are corrections to 
the symmetry factor which arise from the fact that one 
can permute all multiple links joining two vertices 
without changing the diagram, n- is the number of links 
meeting at vertex i of the multi-diagram. 



FIG. 2: (a) A graph consisting of only two nodes 1 and 2. (b) 
The two distinct assignments of the labels 1 and 2 to the linear 
Feynman diagram for I = 4. An example of an inconsistent 
assignment would be for instance 11212 since An = 0. 



III. A FEW EXAMPLES 

Let us illustrate how the method works for spins on a 
graph with N = 2 vertices. We have A12 = A21 = 1, 
All — A22 — and K = 1. As an example we shall 
calculate the contribution of the linear diagram with I — 
4 Hnks from Fig. [TJ Its total contribution is 



1 
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where the first term accounts for the symmetry factor, 
the second one is a product of four propagators, the third 
one is a product of couplings and the fourth one is the 
combinatorial factor P depending on the details of the 
underlying graph's structure encoded in the adjacency 
matrix A. There are only two possible assignments of 
labels 1 and 2 to this diagram, see Fig. [2l All other 
terms vanish and hence P = 2. 

This example is somewhat artificial because the num- 
ber of nodes is small. Now let us suppose that we have a 
graph where pairs of nodes 2j and 2j + 1 are connected, 
so that the only non- vanishing elements of the adjacency 
matrix are A2j,2j+i — ^2j+i,2j = 1 or j = 1, . . . , N/2. It 
is a 1-regular graph. One can easily see that the number 
of decorations P of any Feynman diagram is P = iV or 
P = 0. Indeed, vertices of a given diagram can be alter- 
nately decorated with two consecutive numbers 2j and 
2j + 1. Thus the choice of a single label specifies auto- 
matically the whole decoration. Because this label may 
assume N values, we have P — N. But if the diagram 
has a loop of odd length, one cannot alternately decorate 
vertices along this loop, so in this case P = 0. 

Because P — N, one can see that also F{g) is pro- 
portional to N and hence also In(iVs) ~ N. The pro- 
portionality coefficient /, can be determined in this case 
analytically by a straightforward calculation of the inte- 
gral ©: 



N 

J^(5) = yln 



dx 



'/^erfc^ 



2-2g 



(19) 



with erfc(a;) — 1 — erf(a:). We shall use this explicit result 
to test our method. Using the formula (fT7|) and setting 
K = 1 we find that the first coefficient in the expansion 
of F{g) comes from only one diagram, a line with / = 1 
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from Fig. [11 and reads: 

sPcI^In (V2ttY — = 0.31831... -A^. (20) 

2 V / TT 

The second coefficient /2 is a sum of the previous diagram 
with doubled line, and the one for I — 2: 




On the other hand, we can calculate fis numerically from 
the analytic formula l|19p : 



(22) 



= 0.318315 - 0.1013215^ + 0.0960545^ 
-0.054306/ + 0.055831/ - 0.037248^^ + .... 



We see that both fi and /2 agree perfectly with those 
obtained before. Higher coefficients /; can be calculated 
by performing the sum (flTl) on a computer. We checked 
that all coefficients, up to Z = 11, agree with those from 
Eq. 

In the next section we shall make another cross-check 
by comparing our results with those for the celebrated SK 
model. We shall see that the method produces correct 
values of the coefficients /; as well as of the limiting value 

/* =i^oo(l). 



IV. THE SK MODEL 

The SK model [3| is the spin glass |(T|) on a complete 
graph, where each node is connected to all other nodes. 
Before we apply our procedure of diagrams' summation, 
let us recall what can be calculated for the SK model 
using another methods. The integral |[7]) can be done in 
the thermodynamic limit leading to 

Fooig) = ^ - In (l + erf(gt/V2)) , (23) 
where t is a solution to a saddle-point equation 



t (l + edigt/V2) 



exp 



(24) 



with = Foo(l) = 0.199228.... One can find the coef- 
ficients // by applying Cauchy's differentiation formula 
and by integrating Eq. l(23|) numerically. This gives: 



(25) 



F^{g) = 0.318315 - 0.2026425^ + 0.1474635^ 
-0.1154395^^ + 0.0946265^ - 0.0800585^ + ... 



Let us now calculate /;'s using the method described in 
Sees. II and III. The propagator My for a complete graph 
with TV -I- 1 vertices is = 1/A^ for any pair of i 7^ j. It 
is a -ftT-regular graph with K = N, so as one can see from 
Eq. lfT6|) . each link introduces a suppression factor N~^. 



On the other hand, the combinatorial factor P contains 
a power iV", where v is the number of vertices of the 
diagram. Thus in the thermodynamic Hmit, a diagram 
with I links and v vertices gives a contribution ~ N'"~K 
The exponent v — I is equal to one minus the number 
of closed loops in the diagram. Therefore, in the limit 
N ^ 00 only tree diagrams give non-vanishing contri- 
butions. Our task simplifies therefore to summing only 
tree graphs. Each tree with I links gives the following 
contribution to ff. 



sg 



i+i 



(26) 



where s is its symmetry factor and n.^ is the degree of 
vertex i. 

We performed the summation of all tree diagrams up 
to Z = 11 on a computer and checked that the values 
fi's obtained in this way agree with those obtained from 
Eq. (|23l) . Again, we see that our method gives correct 
coefficients /;. 

We should, however, remember that our goal is not 
only to find the coefficients of the expansion but rather 
/* = Foo{l) which is a sum of infinitely many coeffi- 
cients. If we naively terminate the series at some L: 
Fl{9) = Eti/iff', then e.g. Fn(l) - 0.220701... is 
far away from the true value 0.1992..., because the se- 
ries is slowly convergent. Therefore, we need to find a 
method which allows to read off the limiting value /, 
from first few coefficients. As we shall see below one can 
find a very good estimate of the Hmiting value Foo(l) us- 
ing some general information about the properties of the 
series expansion. Let us first observe that the integral 
([7]) and thus also ^00(5) is convergent only if the matrix 



(27) 



is positive definite. In the Appendix A we show that it 
is so for I5I < 1 and for g = 1, and that G acquires a 
zero mode for 5 = —1, so the integral ^ is divergent for 
5^—1. Prom this we can conclude that the asymptotic 
behavior of the coefficients fi has the following form: 



fi = 



la ■ 



(28) 



with all o; > 0, < a < 1 and limi_^oo ai const. From 
the asymptotic behavior i(28|) we can deduce the following 
formula for i^oo(l): 



Foo(1)-Fl(1) + /l. (^1^ X 
x[C(a,l + L/2)-C(a,l/2 + L/2)], 



(29) 



where C{a,l3) = J^kLoi^ + f^)~" ^ generalized Rie- 
mann Zeta function, and a is estimated from the last 
two coefficients: 



TL-l 



L - 1 



(30) 
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With the help of formula (|29l) we can predict now the 
value of Foo(l) « 0.199226 for the SK model. To esti- 
mate the maximal error we used a method described in 
Appendix B. Our final result /* — 0.199226(5) is in an 
excellent agreement with the analytic result cited above. 



V. RANDOM iC-REGULAR GRAPHS AND 
CAYLEY TREES 

A if-regular graph is a graph with all degrees equal to 
K. We say that a regular graph is random if its adjacency 
matrix A is maximally random under the constraints that 
An — 0, Aij — Aji and J2j — ^ If we 

fix K and let — > oo, random graphs become sparse 
and look locally like Cay ley trees with degree K, because 
the average density of finite-length loops goes to zero in 
this Hmit. Thus for large N, instead of computing the 
coefficients // by averaging them over many if-regular 
random graphs one can calculate them for a single Cayley 
K-tiee. 

The propagator My is simply 1/K if nodes i,j are con- 
nected, and zero otherwise. Unlike for a complete graph, 
the contribution from diagrams with loops cannot be ne- 
glected. Diagrams with multiple connections also do not 
vanish. In order to compute the contribution of a simple 
diagram 7 and its multi-linked versions we need to eval- 
uate the sum lfT7|) . The summation can be performed on 
a computer with only a sHght compHcation as compared 
to the SK model. 

The hardest task in the above is to compute P. Re- 
call that for a 1-regular graph, P was either N or zero. 
Now P is simply equal to the number of different ways 
a given Feynman diagram 7 can be lied down on a K- 
tree, in such a way that any two neighboring vertices of 
the diagram are also neighbors on the tree. One expects 
that the number of such possibilities is of order NK^~^. 
The factor N comes about since one can lie down one 
vertex of the diagram anywhere on the tree. But the 
second has to be located on one of the K neighbors of 
the first one. So each time by finding a position of a 
vertex we should check K neighboring nodes on the tree, 
which gives the factor K^^^. To illustrate this, consider 
again the linear graph for I = 4 from Fig. [TJ The di- 
agram has i; = 5 vertices. The first vertex can be put 
anywhere on the tree, but the next one only on one of 
K neighbors of the first vertex. The same for the next 
one, for which we have again K possibilities, etc., so al- 
together we have P — NK'^. Let us consider now a more 
complicated example of a square graph from Fig. [1] and 
K = 3. Again we have a trivial factor N for choosing the 
position of the first vertex on the 3-tree. Once it is cho- 
sen, we can position remaining vertices only in 15 ways 
as shown in Fig. [H So we have P — 15N which is less 
than P = 3'*A^, because of the constraint coming from a 
closed loop. Clearly, if a Feynman diagram has no loops 
then P = NK'"~^ = NKK For a diagram with loops one 
has to consider constraints on possible decorations. This 




FIG. 3: Example of a Feynman diagram (a square on the top) 
lied down in two different ways on the Cayley tree with K = 3. 
Four links are drawn with four different lines and arrows, in 
order to show how the diagram is put on the graph. When 
one vertex of the diagram (say the upper-left) has fixed index 
i, there are only 6 ways of drawing the diagram as in (a), and 
9 ways as in (b) case. The picture (a) implies that j = I and 
there are 3 ways of choosing j and 2 of k, while for (b) we have 
k — i and 3 possibilities for each index j, I. There is no other 
way of distributing indices j, k, I, so P/N = 2- 3-|-3-3 = 15 
for the square diagram on the 3-tree. 



can be done by enumerating all possible labellings and 
accepting only those which agree with graph's structure. 
It can conveniently be done by a computer program. One 
point must be clarified here - since all P's share the same 
trivial factor TV coming from N possibilities of labeling 
the initial vertex, in the computer program one can just 
fix one vertex of the diagram to have some arbitrary la- 
bel, and consider only decorations consistent with this 
choice. Next, one multiplies the result by A'^, which then 
cancels in the definition of Fooig) that only numbers 
independent on TV remain. 

Using this method, we calculated P's and then the co- 
efficients // up to the given order L (typically L = 11) for 
K = 2,3,4: and 6. The case K — 1 has been analyzed in 
Sec. III. The results are summarized in TableUwhere we 
give the values of /* for K >2 with estimated errors and 
compare them with those obtained by numerical simu- 
lations based on enumeration of all metastable states as 
described in Ref. Q. To save space, the coefficients // 
are not shown and can be found elsewhere 

The agreement with simulations is perfect. For the 
case K = 2 there is also a beautiful analytic result [3| to 
compare with, which gives /* = InA/ir « 0.24156.... As 
we see from the table, it agrees with our result within the 
error bars. We observe also that the rate of convergence 
of the series expansion for Poo (5) grows with K. In other 
words, for smaller K one should go to larger order L 
which is however limited by the fast growth of the number 
of diagrams. This effect is slightly compensated by the 
fact that the complexity of computing the combinatorial 
factor P ~ is smaller for smaller K. 
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VI. REGULAR d-DIMENSIONAL LATTICES 

In this section we shall discuss how to calculate /, for 
d-dimensional regular lattices, namely for d — 2 (square 
lattice) and d — 3 (cubic lattice). The lattices are ir- 
regular graphs with K = 2d, but very special ones. The 
case d = 1, that is a closed chain, gives the same result 
as a random 2-regular graph, because the latter always 
contains at least one long chain, whose contribution to 
Ns in the limit N ^ oo dominates over the contribution 
coming from shorter chains. In the general case, the only 
difference as compared to random regular graphs is that 
now, while calculating P, we shall lie down Feynman 
diagrams on the c?-dimensional lattice. Again we see that 
P is proportional to N since the first vertex can be put 
anywhere, but once it is fixed we can put a neighboring 
one on a neighboring site of the lattice and repeat it 
iteratively for all remaining vertices of the diagram. The 
calculated values of /* are given in Table [J 

VII. LADDERS 

As a further example we consider another particular 
type of if- regular graphs: graphs which we shall call lad- 
ders. A ladder is a graph obtained by stacking above each 
other infinitely many copies of a {d— l)-dimensional cube 
so that corresponding vertices of the copies are aligned 
on a line and the corresponding vertices of consecutive 
copies are joined by a link. A d = 2 ladder is just what 
one usually would call a ladder except that it is infinitely 
long. A d = 3 ladder looks like a bookstand with square- 
shaped shelves. Such ladders are if- regular graphs with 
K = d + 1. It is interesting to compare /* — i^oo(l) 
for ladders with those for random graphs and regular d- 
dimensional lattices. We adopt the usual scheme. The 
only thing which changes is again P, because now we have 
to lie down diagrams on the ladders. The final results are 
presented in Table [B 



VIII. CONCLUSION 

We presented a semi-analytic method to estimate the 
rate /* of exponential growth of the number of metastable 
states Ns of the Ising spin glass on different kinds of 
graphs. We checked that the method reproduces results 
known analytically, which are available for a few par- 
ticular cases. The method is based on a diagrammatic 
representation of the quantity hi{Ns){g) and on the ex- 
act enumeration of all Feynman's diagrams up to a given 
order L. The accuracy of the method improves with 
growing L but already for L close to 10 it yields very 
precise estimates whose uncertainty varies in the range 
of order 10"'' - IQ-^ depending on K (see TableHl. The 
results show that the exponent /, is determined mainly 
by the degree K. This suggests that the number of lo- 
cal minima Ng depends strongly on local properties of 



Graph 


K 


L 


/, , the slop 


e of ln(iV.) 








result 


simulations 


2-reg. graph 


2 


13 


0.2414(2) 


0.242(2) 


3-reg. graph 


3 


12 


0.22484(4) 


0.226(1) 


2D ladder 


3 


12 


0.22568(6) 


0.226(2) 


4-reg. graph 


4 


11 


0.21762(2) 


0.219(1) 


2D lattice 


4 


11 


0.21808(2) 


0.219(2) 


3D ladder 


4 


11 


0.21799(2) 


0.220(4) 


6-reg. graph 


6 


10 


0.21101(2) 


0.211(1) 


3D lattice 


6 


10 


0.21125(1) 




SK model 


oo 


11 


0.199226(5) 


0.199(1) 



TABLE I: Values of /* = i^(l), calculated for various graphs 
with average connectivity K, compared to computer simula- 
tions. L stands for the number of coefficients /; used to esti- 
mate i^oo(l). All simulations were made for A'^ = 10, ... , 24, 
so for relatively small systems. The uncertainty of the last 
digit is given in brackets (). It was estimated as the standard 
error in case of computer simulations, and as the maximum 
error in the way presented in Appendix B for semi-analytical 
calculations. 



the graph and weakly on its global topology. In other 
words, an important information about the complexity 
of the energy landscape of the corresponding spin glass 
is encoded in the short-range properties of the graph. It 
would be interesting to test if this also holds for other 
complex systems. 

The method presented here can be applied to any type 
of graphs. However, the computational complexity of the 
method and the dependence of the accuracy of /* on the 
order L has to be tested case by case. 

In this paper we calculated {Ns) for Gaussian J's. 
Comparing the results to those for binomial J's Q we 
see that (Ns) significantly depends on the distribution 
of J's. It would be quite interesting to investigate the 
dependence on the distribution of J's in a systematic 
way by calculating (Ns) for some other continuous dis- 
tributions of J's. One can try to do this by applying 
the Tanaka-Edwards idea to distributions of the type 
p(J) = / daf{a)e~'^ Another very challenging prob- 
lem is to calculate higher moments (Ns) and eventually 
also {\nNs). 
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Appendix A 

In this appendix we shall prove that the matrix G is 
positive definite for \g\ < 1. We must show that 

X Gx — ^ ^ x-iGijXj , (^■^) 

is strictly positive for any vector x with non-zero length. 
From Eq. (|27|) we have: 

x^Gx^J: . (32) 

Introducing new variables y^: Xi = \/2kiyi we can rewrite 
the right-hand side as: 

2(1 ~g)Y. hyl +gY. A,, (y, + y^f (33) 
- 2(1 +9)Y. hyl + 5 ^ Uj (y,; - y,)^. (34) 

The formula in the upper line shows that the quadratic 
form (|3T1) is non-negative for any x and for < g < 1. It 
is zero only if y = 1 and x = 0. Therefore, F{g) tends 
to a constant for y ^ 1. On the other hand, the second 
equation l(34|) tells us that for —1 < y < it is posi- 
tive definite as well. However, for y = —1 we see that a 
zero mode appears. Because of the zero mode the Gaus- 
sian integration in Eq. ([31 is divergent and consequently 
F{g) ^ oo for y — > —1. These two observations indicate 
that the radius of convergence of F{g) given by Eq. ifTSj) 
is one, and that /; has the form: 

-a;(-l)'r", (35) 

where a € (0, 1] and ai tends to some constant for I oo. 

Appendix B 

The estimation of systematic errors like those involved 
in the calculation of Foo(l) from Eq. l(29|) is not an easy 
task. We believe, that the method described below gives 
an upper bound on the error of Foo{l)- Let us slightly 
adjust the notation for the sake of clarity of the discus- 
sion. The value of Foo{i) from the left-hand side of Eq. 
(|29| depends on L, for which it has been estimated, so it 
is convenient to keep trace of L in the notation. We 
will denote the value of the estimate by F^\l). Of 
course it not the same as Fl{1) which is on the right- 
hand side of Eq. l(29)l and which means just a sum of 



the first L coefficients fi. The method relies on the fol- 
lowing observation. We want to compute the deviation 

Dl ^ F^\l) ~ Fooil) but we do not know the limiting 
value Foo(l). We can however compute slightly modified 

differences d„ = F^\l)-F^\l) for all n = 2, . . . , L- 1, 



IxlO'V 




n 



FIG. 4: Illustration of the procedure of estimating errors. 
We use as an example the data for the SK model, for L = 
11 known coefficients fi. We plot {n,d„) (squares) for n = 
2, . . . , 10 in a log-log plot, and then fit a straight line. The 
line is shifted to be above all the data points. We extrapolate 
the line to obtain the value for n = L = 11 (circle) and take 
this as an uncertainty of -Foo(l). 



and say that Z?„ ~ d„ for L being sufficiently large. Now 
we can plot c?„ versus n and extrapolate it to n = L to 
obtain an estimator of c?l which in turns estimates Z?l 
giving the error. Because dn falls with n as a power of 
n or faster, we can estimate dh from above by fitting 
a straight line to the points (lnc?„,lnn) and shifting it 
so that all points lie below it. By extrapolating it to 
n = L we get the value dh and use it to estimate the 
upper bound for the deviation between fI^^^ (\) and the 
true value Foo(l)- This procedure is illustrated in Fig. [4] 
for the case of SK model from Sec. III. How rehable is 
this method? We checked that intervals F(\) ± Ai^(l) 
obtained from n = 2, 3, 4, . . . , L — 1 first coefficients fi 
always include the value for n ^ L for all graphs dis- 
cussed in this paper. We checked also that reducing the 
error, say, by a factor two would result in many situa- 
tions for which Fl{1) would lie outside the error bars. 
This means that the method does not overestimate the 
error too much. 
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